Decoherence-free dynamics of quantum discord for two correlated qubits in 

finite-temperature reservoirs 
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We investigate decoherence-free evolution (DFE) of quantum discord (QD) for two initially- 
correlated qubits in two finite-temperature reservoirs using an exactly solvable model. We study 
QD dynamics of the two qubits coupled to two independent Ohmic reservoirs when the two qubits 
are initially prepared in X-type quantum states. It is found that reservoir temperature significantly 
affects the DFE dynamics. We show that it is possible to control the DFE and to prolong the DFE 
time by choosing suitable parameters of the two-qubit system and reservoirs. 
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I. INTRODUCTION 



It is well known that quantum entanglement [!H7|; is 
a distinctive quantum feature of quantum correlations, 
but not the only one. In order to capture nonclassical 
correlations, several measures of correlations have been 
proposed in the literature fc^-flij . Among them, quantum 
discord (QD) has received considerable attention [15l — 
H3] , which quantifies the quantumness of correlations be- 
tween two partitions in a composite state [2l[. The QD, 
which measures general nonclassical correlations includ- 
ing entanglement, is defined as the mismatch between two 
quantum analogues of classically equivalent expression of 
the mutual information. For pure entangled states, all 
nonclassical correlations which can be characterized by 
the QD arc identified as entanglement. For mixed states, 
the entanglement can not totally describe nonclassical 
correlation. Even for some separable states, although 
there is no entanglement between these two parts, the 
QD is nonzero, which indicates the presence of nonclassi- 
cal correlations. And compared to classical counterparts, 
such separable states can also improve performance in 
some computational tasks [22|, [23[. In fact, Nonclassical 
correlation described by the QD can be considered as a 
more universal quantum resource than quantum entan- 
glement in some sense, and the QD offers new prospects 
for quantum information processing. 

Any realistic quantum system will inevitably inter- 
act with the surrounding environment, which causes the 
rapid destruction of crucial quantum properties. There- 
fore, besides the characterization and quantification of 
correlations, an interesting and crucial issue is the be- 
havior of correlations under decoherence. The QD dy- 
namics have been widely studied by using models where 
qubits interacts in zero-temperature reservoirs [2(| [U- 
[301 ]. Moreover, QD is more robust than the entan- 



glement against decoherence [24|, 1311-33] . Also, unlike 
the entanglement, which exhibits sudden death (33 - [39| . 
the QD displays the sudden-transition phenomenon from 
a decoherence-free-evolution (DFE) to a decoherence- 
evolution regime in the dynamic evolution p6l - [28} . and 
the inital QD remains unchanged in the DFE regime. 
The This sudden-change phenomenon have been demon- 
strated in recent experiments [29, 30]. It is of significant 
interest to prolong the DFE time of the QD in quantum 
information science since quantum information process- 
ing favors the long DFE time of the QD in the dynamic 
evolution. In this paper, we study the possibility of pro- 
longing the DFE time of the QD by investigating the 
DFE dynamics of initially-correlated two qubits in two 
finite-temperature reservoirs. We shall show that the 
DFE time of the QD can be controlled by choosing initial- 
state parameters of the two qubits, reservoir parameters, 
and qubit- reservoir interactions. 

This paper is organized as follows: In Sec. [Til we 
present our physical model and study its solution. In 
Sec. IIII1 we study the DFE dynamics of the QD to indi- 
cate the possibility of prolong the DFE time of the QD by 
investigating QD dynamical behaviors of the two qubits 
in finite-temperature reservoirs for certain initially pre- 
pared A-type states. The effects of the temperature, the 
initial-state parameters, the system-reservoir coupling on 
QD are studied in Sec. lIVI Finally, we conclude this work 
in the last section. 



II. 



THE MODEL AND ITS SOLUTION 



'Author to whom any correspondence should be addressed; Elec- 
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The total system we consider is shown in Fig. [TJ in- 
cluding a pair of noninteracting qubits A and B in two 
independent heat baths, respectively. The energy separa- 
tion between the excited state \e) and ground state \g) of 
qubits is denoted by fit (i = A,B). Atoms A and B are 
coupled individually to its thermal bath with tempera- 
ture T a and respectively. Here, each bath is modelled 
as an infinite number of harmonic oscillators with fre- 
quencies uj a k and ujf,k, which couple to the relevant elec- 
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FIG. 1: (Color online) Schematic of a pair of noninteracting 
qubits A and B with energy separations Qa and Q,b under the 
influence of two independent environments with temperature 
T a and T b . 



tronic degrees of freedom of i-th qubit via the coupling 
constants gi k . The Hamiltonian of the total system in- 
cluding the two qubits and the environment is composed 
of three parts 



H = Hs + Hb + Hse j 



(1) 



where Hamiltonians of the system and baths are given 
by 



H - QA ~z J- Qb ^ 
H S - —°A + ~ a Bi 



Hp 
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^bkb[b k 



(2) 



(3) 



where of is the standard diagonal Pauli matrix, and 
a k (b k ) and a k (bk) are the creation and annihilation 
operators for a oscillator in fcthe mode, obeying the 
bosonic commutation relation. The system-bath inter- 
action Hamiltonian is given by 



Hse — ft {u z A 9Akak + cr%gBkbk 



■h.c. 



(4) 



Ui(t) = exp 




XkCxk W-h-c] 



(6) 



where £ X k (*) 
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) . Here, X = a(b) when 



i = A(B). The analytical solution of this model allows us 
to study behaviors of these correlations under the action 
of decoherence without any approximations. 

In order to show how decoherence affects the corre- 
lations in this two-qubit composite system, we assume 
that the initial state of the two baths is a thermal state 
denoted by the following density operator 



Mo)=n 



i 



-0xHuixk\ p-pxhuxkXtXk 



(7) 



X.k 



where fix = {ksTx) , and the initial state of the two 
qubits is denoted by the density operator ps (0), the total 
system is assumed as a product state of these two initial 
states p(0) = ps (0) ® pe (0). By tracing out over the 
state of the environment, we can obtain the quantum 
state of the two qubits at any time 



ps (t) = Tr B [tft (t)p(0)U(t)] 



(8) 



In the Hilbert space spanned by the two-qubit product 
state basis {\gg) , \ge) , \eg) , |ee}}, the density operator 
ps (t) loses its off-diagonal terms as a result of the in- 
teraction with the environment E. The elements of the 
density operator ps (i) can be exactly obtained as 

P(U,lB)(jA,jB)( t ) = P(Ia,Ib)(ja,Jb)( ) 

xexp[(S lA , jA -l)T A (i)] 
xexp[(5 lB j B -l)r B {t)], (9) 

where l,j = e,g and is the Kronecker delta func- 
tion. Here, we note that no approximation is employed. 
Equation shows that diagonal terms of the reduced 
density matrix ps (t) remain in the initial value, however, 
the off-diagonal terms vary with the time evolution. The 
two time-dependent decohering factors arc defined by 



The bilinear interaction between the system and the 
bath in Hsb indicates that the atomic inverse operator 
erf commutates with Hamiltonian H, i.e. [a?,H] = 0. In 
the rotating frame with respect to Hamiltonian Hs+Hb, 
the interaction Hamiltonian reads 

H (*) = J2 h {9Ake- l ^ kt <J A a k + g B ke- l ^ kt n B b k + h . c . ) , 

k 

(5) 

Equation |5]) indicates that the state of the environment 
will be sensitive to the values of erf. The commutation 
[Hs,Hse] = suggests that no energy exchanges be- 
tween qubit and bath, i.e., energy of the system S is 
conservative [4(J IH|. Therefore, the model describes a 
purely decohering mechanism. The evolution operator 
generated by the effective Hamiltonian ([5]) is given by 
U (t) = U a (t) U b (t) with 



X,k 



\Zx,k(t)\' 



coth 



( h(3 x u X k \ 



(10) 



which are two real- value functions completely to charac- 
terize the two-qubit dynamics in the decohering process. 

Since the states in the reservoir are very dense (con- 
tinuum), we can take the continuum limit to convert the 
summation over k into an integral with respect to u)x- 
Then the properties of environment are described by the 
reservoir spectral density 



J (tJ X ) = ^S(u>x - wx,k) \9x,k[ 

X,k 



(11) 



With these replacements, the two decohering factors in 
Eq. (10) become 



r, (*) 



x ^ x iM coth (M^L } (12) 
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III. DFE DYNAMICS OF THE QD FOR TWO 
INITIALLY CORRELATED QUBITS 

In this section we study DFE dynamics of the QD for 
the two initially correlated qubits in the two independent 
finite-temperature reservoirs by investigating the QD dy- 
namics under decoherence. The total correlations be- 
tween two qubits A and B described by a bipartite quan- 
tum state pab — Ps are generally measured by quantum 
mutual information 



l{ps) = S (p A ) + S (p B ) - S (p s ) 



(13) 



where S (p) = — Tr(plogp) is the von Neumann entropy 
of density matrix p, pa (pb) is the reduced density oper- 
ators for subsystem A(B). Quantum mutual information 
contains quantum correlation V (ps) and classical corre- 
lation C (ps)@, H3, EH- Quantum correlation V (ps) can 
be quantified by the QD [lj, [21j , which is defined by the 
discrepancy between quantum versions of two classically 
equivalent expressions for mutual information 



V(p s )=X(ps)-C(ps). 



(14) 



where the classical correlation C (ps) is defined as the 
maximum information about one subsystem pi, which 
depends on the type of measurement performed on the 
other subsystem. For a local projective measurement life 
performed on the subsystem B with a given outcome fc, 
we denote 



Pk = Tr AB [(I A ®n k )p s (I A ®IL k )}. 



(15) 



as the probability, where I A is the identity operator for 
the subsystem A. Then the classical correlation reads 



C( P s 



max 
{n fc } 



S(pj 



(fe)x 
A i 



S(p A ) - min 
{n fc } 



^2pkS(p 

k 



(16) 



L k 



where the maximum is taken over the complete set of 
orthogonal projectors {11^} and 



Here, the coefficients Cj with real constants satisfy the 
condition that ps(0) is positive and normalized, and Iab 
is the identity operator of the two qubits. Obviously, 
such states ps (0) are general enough to include states 
such as the Werner states and the Bell states. 

Using the results obtained in Eq.©, the density ma- 
trix of two qubits at time t has the following analytical 
expression 



(19) 



where a = a(t)e^ nA+QB ^ and 7 = ^ e i{n B -U A )t withj 

a(t) = (c 1 -c 2 )D A (t)D B (t), (20a) 
7(t) = (c 1 +c 2 )D A (t)D B (t), (20b) 

where Di(t) = exp [— (t)] are decohering functions. 
Since the reduced density matrix pi (t) — i^/2 of sub- 
system i = A,B, and the eigenvalues of the two-qubit 
density operator ps (t) can be exactly calculated from 
Eq. (fT9|) with the following expressions 
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l + c 3 



Ai. 



(1 + c 3 T a) 



^3,4 = T (1 - C 3 =F 7) , 

quantum mutual information can be written as 



l(p s ) =2 + ^AilogAi. 



(21a) 
(21b) 

(22) 



Now, we turn to calculating the classical correlation 
C(ps) defined in Eq. P^)l . In the Hilbert space spanned 
by {|e) , \g)}, any two orthogonal states |0) and |1) can 
be represented as a unitary vector on the Bloch sphere 



|0) =cos%) + e^sin0|e), (23a) 
|1) =e- # sin%) -cos0|e), (23b) 



pf - -Tr s [(Ia 
Pk 



'Hk)p(iA ®n fc )] 



(17) 



is the reduced density matrix of subsystem A after ob- 
taining the measurement outcome k, which is normalized. 

Let us first find the analytical expression for the den- 
sity operator of two qubits at time t. We assume the 
two qubits are initial in a class of states with maximally 
mixed marginals, which is described by the X-structured 
density operator 



with < 6 < tt/2 and < < 2n. Therefore, for a local 
measurement performed on the subsystem B, a complete 
set of orthogonal projectors contains two elements IIo = 
0) (0| and III = |1) (1| with the probability Po = Pi — 
1/2. After the two project measurements the reduced 
density operator of subsystem A with an outcome k reads 



Pa 



2(1 -c 3 cos(26>)) (-1) esin(20) 
(-l) fe e*sin(26>) 2(1 + c 3 cos(20)) 



(24) 



Ps (0) 



Iab 
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E 



Cj^A 



(18) 



where we have introduced the following parameter 

e = ae i[{^A+n B )t-<t>] + rye i[(n A -n B )t+,p]_ ^5) 
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Making use of Eq. (24), it is straightforward to calcu- 
late the classical correlation with the following expression 

M 



C (/*(<)) = £ 



1 + (-!)"*(*) 



log 2 [1 + (-lTx(t)} , 

(26) 



where 



X (t) = max[|c 3 |,(|a(t)| + | 7 (i)|)/2] 



(27) 



which depends on the relation between the coefficients Cj 
of the initial state in Eq. (|18|) and the dynamic parameters 
defined by Eq. (22). From Eq. (27) we can see that: 
(i) If |c 3 | > (M + | 7 |)/2, we have X {t) = |c 3 |; (ii) If 
|c 3 < (|a| + |7|)/2, we have x(t) = (\a\ + |-y |) /2. 
Then the QD between the two qubits can be written 

as 



where we have introduced £l±(t) = 1 ± DA(t)Ds(t). 

Subtracting the classical correlation from the quantum 
mutual information I(ps), we find that the QD is given 
by 

1 i 

n=0 

which is independent of time. Hence, in the regime of 
t < t p the QD dynamics is decoherence-free, and the 
initial QD can be preserved. 

In the regime of the QD decay, i.e., t > t p , since 
DA{t)DB{t) < |c 3 |, we have x(i) = |c 3 |. The classical 
correlation is a constant 



C(ps(t)) = E ^[1 + (-l)"c 3 ] log 2 [l + (-l)"c 3 ], (34) 

?i=0 



V(p s (t)) = 2 + E A « lo 82 ^ - C(ps(t))- (28) th en the QD reads 



i=l 

In order to clearly understand the QD dynamics of two 
qubits in two independent reservoirs, we consider the case 
of ci = 1, C2 = — c 3 , and |c 3 | < I. In this case, from Eq. 
(f2~T| we can obtain 

X (t)=max[\c 3 \,D A (t)D B (t)], (29) 

where the two decohring functions are given by 



Aw=(i+^ 2 r*n 



n=l 



2 j-2 



1 



UJtt 



-I ~Vi 



(I + h(3 i u> ic n)' 1 



( 3 °) 

where i = (A, B) and we assume the two reservoirs have 
the Ohmic spectral density [43[ given by 



J{wi) = ^Wjexp 



(31) 



(35) 

which is a decaying function of time. 

The above discussion shows that the QD exhibits a 
sudden change at the critic time t p . The QD evolution 
dynamics displays different behaviors before and after the 
critic time. The initial QD can be preserved in the regime 
of t < t p while the QD decays in the time evolution 
regime of t > t p . In what follows we discuss the depen- 
dence of the DFE time of the QD upon the initial-state 
parameters, the system-reservoir coupling, and the reser- 
voir temperature for the case of two equal-temperature 
reservoirs and the case of two unequal-temperature reser- 
voirs, respectively. 



where rji > is the system-reservoir coupling constant, 
and LOi C is the high-frequency cut-off. 

Taking into account the decohring functions DA{t) 
and Dsit) decay in the time evolution, making use of 
Eqs. (21), (26), (28), and (29), we can find the critic 
time denoted by t p which is given by the condition 
DA(tp)DB(t p ) = |c 3 |, and we can see that before the 
critic time, i.e., t < t p , the initial QD can remain un- 
changed in the time evolution while the QD decays after 
the critic time, i.e., t > t p . Therefore, the time t p is 
the DFE time of the QD. In what follows we calculate 
explicitly analytical expression of the QD in the whole 
time evolution to indicate how to control the DFE time 
of the QD. 

In the DFE regime of the QD, i.e., t < t p , since 
D A (tp)D B (t p ) > \c 3 \, we have x (f) = D A (t)D B (t) = 
DaDb- Consequently, the classical correlation reads 
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FIG. 2: (Color online) 3D diagram of the quantum discord as 
a function of /3 and u) c t in the case of two equal-temperature 
reservoirs under the setting r\ — 0.2, ci = 1, and C2 = — o A = 
0.4. 



(32) 
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A. The case of two equal-temperature reservoirs 

In this subsection, we assume the two reservoirs are 
identical, i.e., uj Ac = w Bc = uj c ,r/ A = i] B = 77, j3 A = /3 B = 
P. Then we have D A (t) = D B (t) = D(t). 



the DFE time of the initial QD t p increases and the QD 
decay after the time t p becomes slower when the system- 
reservoir coupling constant 77 decreases. Hence, we can 
conclude that the smaller the qubit-reservoir coupling 77, 
the more robust the QD against decoherence. 
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FIG. 3: (Color online) Plot of the quantum discord vs the 
scaled time uu c t in the case of two equal-temperature reservoirs 
with the qubit-reservoir coupling 77 — 0.2 for dash green line, 
r\ — 0.6 for solid blue line, 77 — 0.9 for dash-dot red line. Other 
parameters are taken as ci = 1, c% = — C3, and C3 = 0.4, and 
= 5. 



FIG. 4: (Color online) Plot of the quantum discord vs the 
scaled time ui c t in the case of two equal-temperature reservoirs 
with the initial-state parameter | C3 1 — 0.2 for dash green line, 
I C3 1 = 0.4 for solid blue line, and |ca| = 0.8 for dash-dot red 
line. Other parameters are taken as 77 = 0.2, ci = 1, C2 = — C3, 
and /3 = 5. 



At zero temperature, from Eq. (30) we can explicitly 
obtain the decohering function with the form D(t) — 
(1 + uj 2 t 2 )~ v / 2 , which allows us to achieve the expression 
of the QD preservation time 

^ = , (ae) 

which indicates that the DFE time depends on the initial- 
state parameter C3, the qubit-bath coupling 77, and the 
bath parameter uj c . The DFE time increases with the 
decrease of the initial-state parameter | C3 1 or/and the de- 
crease of the qubit-bath coupling 77. Moreover, for a set of 
given parameters (03,77), the lower the cut-off frequency 
of the reservoir is, the longer the preservable time will 
be. 

At finite temperatures, we can numerically study the 
influence of the initial-state parameter, the qubit-bath 
coupling, the bath parameter, and the bath temperature 
upon the QD preservation time. In Fig. [21 we have 
plotted the QD as a function of (3 and Lu c t. From Fig. 2 
we can see that a sudden change of the QD occurs at the 
critic time of t = t p , the initial QD can be preserved in 
the time evolution regime of < t < t p . However, The 
QD decays in the time evolution regime of t > t p . We 
can also see that the decrease of the bath temperature 
prolongs the DFE time of the initial QD and slows the 
decay rate of the QD. 

Fig. [3] shows the effect of the qubit-reservoir coupling 
77 on the QD dynamics. From Fig. 3 we can see that 



The effect of the initial-state parameter | C3 1 on the QD 
dynamics is shown in Fig. 0] where we have plotted the 
time evolution of the QD for (3 = 5. From Fig. 4 we can 
see that one can make the DFE time of the initial QD 
t p longer by decreasing the magnitude of the initial-state 
parameter | C3 1 . 



B. The case of two unequal-temperature reservoirs 

In this subsection, we the influence of the temperature 
difference the QD dynamics when the two reservoirs have 
different temperature. We suppose that uj Ac — uj Bc — 
^>c,tia = f]B = T), and the temperature of two reservoirs 
satisfies the relation T a — nTf,, i.e., /3& = n/3 a where k is a 
parameter to denote the temperature difference between 
the two reservoirs. In case of two unequal-temperature 
reservoirs, the QD dynamics is described analytically 
Eqs. (33) and (35). In Fig. 5, we have plotted the 
dynamic evolution of the QD with respect to f3 a for dif- 
ferent values of the temperature-difference parameter k 
when the two qubits are initially prepared in the A-type 
state with the state parameters c\ = 1, Ci = — C3, and 
c 3 = 0.4. 

From Fig. 5 we can see that the temperature difference 
between the two reservoirs does not change the maximal 
value of the QD but does change the DFE time of the 
QD in the dynamic evolution. The maximal value of the 
QD is stable for different temperature-difference param- 
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FIG. 5: (Color online) 3D diagram of the quantum discord 
as a function of f3 a and uj c t in the case of two unequal- 
temperature reservoirs under the setting: r\ — 0.12, ci = 1 
and C2 = — C3 = 0.4 for different value of n, (a) n = 0.2, (b) 
K = 1, (c) K = 5. 



ronments through investigating the QD dynamics of the 
two qubits coupled to independent reservoirs. We have 
chosen the initial state of the two qubits to be X-type 
states which is known to exhibit quantum correlations 
described by the QD. We have shown that reservoir tem- 
perature significantly affects the DFE time of the QD. 
We have also indicated that it is possible to control the 
DFE dynamics of the two qubits and to prolong the DFE 
time by choosing suitable parameters of the two-qubit 
system and their environments. These results shed new 
light on the QD control which would be a new direction 
in quantum information processing. 

For X-type initial states with the state parameters 
ci = 1, C2 = — C3 and | C3 1 < 1, we have investigated 
the QD dynamics of the two qubits in detail when the 
two independent reservoirs are Ohmic. In the case of 
the two zero-temperature reservoirs, we have shown that 
the DFE time increases with the decrease of the initial- 
state parameter | C3 1 or/and the decrease of the qubit- 
bath coupling rj. Moreover, for a set of given parameters 
(03,77), the lower the cut-off frequency of the reservoirs, 
the longer the DFE time is. In the case of the two nonzero 
equal-temperature reservoirs, we have indicated that the 
decrease of the bath temperature can prolong the DFE 
time of the initial QD and slows the QD decay rate. In 
the case of the two unequal-temperature reservoirs, it is 
found that the temperature difference between the two 
reservoirs does not change the maximal value of the QD 
but does change the DFE time in the dynamic evolu- 
tion. The larger the temperature difference parameter, 
the longer the DFE time t p . In this sense, we conclude 
that the DFE time can be prolonged under certain condi- 
tions through the increase of the temperature difference 
between the two reservoirs. 
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